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In conducting ferromagnetic materials, a moving domain wall induces eddy currents in the sample 
which give rise to an effective retarding pressure on the domain wall. We show here that the pressure 
... is not just proportional to the instantaneous velocity of the wall, as often assumed in domain wall 

' models, but depends on the history of the motion. We calculate the retarding pressure by solving 

, the Maxwell equations for the field generated by the eddy currents, and show how its effect can 

' be accounted for by associating a negative effective mass to the magnetic wall. We analyze the 

(N ■ dependence of this effect on the sample geometry and discuss the implications for Barkhausen noise 

(~| , measurements. 

PACS numbers; 75:60.Ch, 75:60.Ej 

I. INTRODUCTION 



Soft magnetic materials subject to a slowly changing external magnetic field, respond with a jerky motion of the 
domain walls, known as Barkhausen noise Many universal properties of this noise have been identified, and are 

^ ' correctly reproduced by theoretical models, confirming that, in analogy to critical phenomena, the statistical features 
of the signal only depend on general properties of the physical mechanism that govern the magnetization reversal 
process, while they are independent from the microscopic details of the particular system [3, 0, Q. 

The average shape of the pulse, being a non-scalar quantity, has been successfully identified as a powerful tool to 
characterize the universal properties of crackling systems: pulses of different duration, once properly rescaled, are 
' expected to collapse onto a universal function_[6| . In the case of Barkhausen noise, pulses from experimental data do 
approximately collapse on the same curve 0, [J, however, this curve shows a clear leftward asymmetry, while models 
Q that very accurately reproduces most of the other universal quantities predict a symmetric shape Q . 

, As shown in Ref. the asymmetry in Barkhausen pulses is due to the non-instantaneous response of the 

eddy current field to the domain wall displacement. Barkhausen noise models usually assume this response to be 
I [ instantaneous, and thus do not capture this asymmetry. Since eddy currents take a finite time to set up and also 
^ 1 they persist for a finite time after the corresponding wall displacement, the pressure on the moving domain wall is 
fSl ' not strictly proportional to the instantaneous velocity of the wall, but depends on the history of the motion. The 
' delay has a characteristic timescale, and therefore its effect is more evident on avalanches of comparable duration, 
' and disappears on very large ones, where the separation of timescales is such that the response of the field can be 
, assumed to be instantaneous, and strict universality is recovered. The first order correction to the instantaneous 
' response approximation can be accounted for by associating a negative effective mass to the wall in the equation of 
O ' motion fvq. 

In this paper we report a detailed calculation of the retarding pressure starting from Maxwell equations, and obtain 
' the negative effective mass as a first order correction to the quasi-static approximation. Our approach is similar to 
^ . the one of Bishopfioj. The resulting non-local damping was previously employed in a domain wall dynamics model 
' in Ref. where its effect on the Barkhausen pulse shape was studied and compared with experimental data. Here, 
, we analyze the role of sample geometry on the eddy current propagation and provide expression for the damping term 
d ' and the effective mass as a function of the sample aspect ratio. 

The paper is organized as follows. In Sec. II we solve the Maxwell equations for the eddy currents in a conducting 
sample with a moving domain wall. In Sec. Ill we compute the resulting pressure on the domain wall. In Sec. IV we 
derive the first order correction to the pressure. In Sec. V we discuss the role of the sample geometry and in Sec. VI 
we conclude. Finally, two appendix report the details of some series summation used in the manuscript. 



II. THE EDDY CURRENT FIELD FROM THE MAXWELL EQUATIONS 

Consider a sample with dimension x £ [—a/2, a/2], y G [—6/2,6/2], and infinite in the z direction, divided in two 
magnetic domains by a rigid domain wall on the yz plane, moving from position a; = 0, as in Fig[TJ The displacement 
of a magnetic wall in a conducting medium induces a flow of eddy currents that generates a magnetic field, which, in 
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the geometry indicated in Fig. [Tl is parallel to the z axis: 

H = Hix,y,t)z. (1) 

The magnetic field obeys the Maxwell equation 

V'^H afj.dtH , (2) 

where cr and fi are the electric conductivity and the magnetic permeability of the medium. In Eq. [5] the displacement 
currents are neglected with respect to the ohmic currents. This equation is usually solved in the quasi-static approx- 
imation, where ^ is negligible within domains, and the equation reduces to V^iJ — 0. Eq. [2] is a diffusion equation: 
there is a finite time delay between the wall displacement, and the establishment of the eddy currents. Given that 
the typical timescale for diffusion is proportional to /i, the quasi-static approximation /i = corresponds to assuming 
an instantaneous response of the field. 
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FIG. 1: Horizontal arrows indicate the directions of the magnetization in the two domains. The black vertical arrow indicates 
the direction of motion of the wall. 

To take into account dynamical effects from eddy currents in the domain wall propagation, we need to solve Eq. [21 
with the appropriate boundary condition _ff = on the sample surface. The discontinuity across the wall is regulated 
by the Faraday condition dxH{0'^ ,y, t) — dxH{0^ , y, t) = 2alv{t), where v{t) is the velocity of the wall. 

Let us expand H in its Fourier components 



H{x,y,t) 



/27r 



dujF{x,y,uj)e'" 



Each component has to satisfy 



VF(x,?/,w) = r'^F{x,y, 



(3) 



(4) 



with = iufia. The magnetic field has to be zero on the sample boundary, which implies the conditions 
F{x,ztb/2,uj) — F{±a/2,y,uj) = 0. Moreover, the Faraday condition around the wall imposes dxF{Q'^ ,y,LLj) — 
dxF{0~ ,y,uj) = 2aIv{Lo), where / is the saturation magnetization, and v is the Fourier transform of the velocity of 
the wall v{t). The solution is of the form 



= ^ A„(x, uj) cos(A„y) , 



n=0 



where An satisfies 



with Kl = 



dlAn{x,Ul) = kl^An(x,Uj) , 



(5) 



(6) 



^, to be solved separately for x > and a; < 0. The condition on the sample boundary in the y 
direction implies cos(±5/2A„) = which fixes A„ = (2n + l)7r/5. The condition on the other boundary in satisfied by 
choosing 



An{x,Lu) = C„(t^)sinh(A„(| x \ -a/2)). 



(7) 



so that 



F{x,y,u}) = ^ C„(cj) sinh(A„(w)(| x \ -a/2)) cos(A„y) . 



(8) 



n=0 
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The functions C„(w) are fixed by the Faraday condition 

oo 

d,F{0+,y,oj) = ~d,F{0-,y,uj) = ^ C„(c^)A„(c^) cosh(A„(u;)a/2) cos(A„y) = <jIv{uj) . (9) 



ri=0 



Multiplying by cos(Amy), integrating in [—6/2,5/2], and using the orthogonality relations /^(/2 cos(A„?;) cos(Amy) 
Sn,nib/2, and j'^i%dy cos{Xmy) = (-1)™2/A™ we get 



C„(c^) = (-1)"^^^^^ TT^^^^vi^) , (10) 

^ ^ ^ ^ 6 A„A„(w)cosh(A„(w)a/2) ^ ^ ' 



so that finally 



Fi.,yM = g(_i)nl^ tanh(A W2) 

III. EDDY CURRENT PRESSURE ON THE WALL 

Once the solution of Maxwell equation with the appropriate boundary condition is given, the average eddy current 
pressure on the wall is obtained by integrating the magnetic field over y at the wall position a; = 0: 

21 f''/^ 2T 1 1"°° f^/^ 

P{t)^- dyH,{Q,y,t)^-^ dioe''^' dyF^{Q,y) , (12) 

J-b/2 V 27r J-oa 

or, in terms of the Fourier transform: 

P{iu) = -v{oj)f{uj) , (13) 

with 



it. .\ _ IP-' tanh(A„(^Ja/2) 

In real space the pressure at time t is given by a convolution of velocities of the wall at all times prior to t with the 
response function /: 

P{t) = ^j^ dsv{t-s)f{s). (15) 
To simplify the uj dependence in Eq. [141 we use the relation 

oo ^ 

tanh(.) = 8z g ,2(2fc + i)2 + 4,2 ' (16) 

with z — A„(cj)a/2, which gives 

^ 64a/V 1 

■^^'^^^ 52 ^^^A2(62A2+„2A2(^))- ( ) 

Replacing A„ and A„ with their expressions, we get 



54/2 00 -[^ 

/(w) = — p5 ^ > ^5 7715 ^ , (18) 
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where ^ indicates a summation over odd numbers only, and uia = t^^ = n"^ / a iia^ , uJb = ^ — ■K^/a^b'^. Given that 
l/{ujQ + iuj) has ^/2Tr exp{—uJot)9{t) as inverse Fourier transform (where 9{t) is the Heaviside theta function), it is easy 
to obtain from Eq. [15] the expression for / in real space by anti-transforming term to term in the double sum to get 

n,k—l 

where ujk.n = ~ ^^'^a + n^ujb- The response function results to be the sum of simple exponential relaxations, with 

different relaxation times. The largest and therefore most relevant relaxation time is tq^q = ^ + ^) ^ . 

The form[T9]of the response function gives through the convolution [12] the explicit solution for the retarded pressure 
on the wall. 



IV. FIRST ORDER CORRECTION: DAMPING COEFFICIENT AND NEGATIVE EFFECTIVE MASS 

Given the full solution of the response function, wc now want to calculate the first order correction to the quasi- 
static approximation fj, — in the retarded pressure. In order to do this, let us replace the expression [19] in the 
convolution [12] and exchange the sum with the integral: 



n.k—\ 

where the Heaviside function was eliminated by restricting the domain of integration. For small relaxation times 
(note that rfe^„ < ro,o for every fc and n) the exponential functions in the integral decay very fast around s = 0, so 
that the velocities that sensibly contribute to the convolutions are only those at time very close to t. We may thus 
expand v(t — s) ~ v{£) — sv'{t) around s = t, and perform the integrals. The term proportional to v{t) corresponds 
the usual instantaneous contribution, while the term proportional to v'(t) gives rise to the first order correction to 
the quasi-static approximation: 



Replacing the frequencies with their expression 

P{t) c ^i^^ ivm,ib/a) - v'{t)nJ:2{b/a)) , (22) 
where Si (a) = J2n,k WWlk^^ ^2(0;) = J2n,k ij^T+h^^ ■ Eq. [22] allows to identify a damping coefficient 

^ 64/2cr&2 



and an effective mass 



4 



-Si(V«), (23) 



64/20-2 „54 

M-^ r Mb/a), (24) 

which turns out to be negative. Moreover, a characteristic time r of order /i can be identified as the ratio between 
mass and damping: 

r=m/T^nmfl. (25) 

A similar calculation to get the damping coefficient and the effective mass may be carried out in frequency space. 
From Eq. 1181 that gives the Fourier transform of the full response function, we can separate the real and the imaginary 
part to get 



64/ k (Va + n iOb — 



/('^) = 7I2T:2 ^2^., //,.2.. , ^2..,^2 , ..2^ ' (26) 
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which allows, by writing f{uj) = r(w) — icuM, to formally identify a frequency dependent damping coefficient 



and effective mass 



_ UPa^ ^ 1 



Expanding to the first order in fj. 



/H~ ^^{h/a)~iu J^S2(6/a), (29) 



we recover F and M as in Eqs. I23I24I 



V. DEPENDENCE OF DAMPING AND MASS ON THE GEOMETRY OF THE SAMPLE 

Eqs. I24I23I25] give the expression correct to the first order in /i of damping coefficient, effective mass and charac- 
teristic time in a general geometry, defined by the parameters a and b. The series Si and S2 can be summed up for 
particular geometries. The calculations are reported in the appendixes. 

For a slab with b 

j:,{b/a ^ 0) ^ , Mb/a ^ 0) ^^As , (30) 

where A„ = Efc^"": with A3 ^ 1.05179... and A5 = 1.00452..., so that 

r2 ^^^3;, AT r2 2 8A5 3 A5 2 /q1^ 

1=7 cr — —b , M = —I a LL — r-6 , r — 'Jn-—^——b . (31) 
TT'' TT° 27r^A3 



For a slab with a <^ b 



so that 



For a square rod with a — b 



so that 



I]i(6/a-^c5o)~— — , E2(fo/a^oo)~— — , (32) 



F = Paa , M = -Pa^u— , t = an— . (33) 
^12 ' 12 ^ ' 



Si(l) = _ , S2(l) ~ 0.264 ± 0.001 , (34) 



F = /V^ , M ~ -/VV 0.0175 b^ , Tc^afi 0.0351 5^ . (35) 







a < & 


a = b 


T/{Pax) 


0.054 


1 


0.5 




0.026 


0.083 


0.0175 




0.048 


0.083 


0.035 



TABLE I: Values of F, M, and r for varies geometries, x denotes the smaller between a and b. 
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The results are summarized in Tab. I. The dependence on a (b) correctly disappears in the limit a — > oo (& — > oo). 
The physical quantities F, M , and r all increase with the overall sample size. For example, for a sample of a given 
thickness a, both damping coefBcient and effective mass increase with b and respectively. However, as soon as 
b becomes larger than a, they saturate to a value proportional to a and a'^ respectively. Since the mass increases 
faster than the damping with b, the characteristic time also increases with b, and saturates to a value proportional 
to a^. A similar behavior is observed by varying a at fixed 6, although the role of the two dimensions transverse and 
parallel to the wall is not symmetric. Essentially the dependence of F, M, and r on the geometry of the sample is 
dominated by the smaller between a and b. Thus the relevance of the eddy current dynamic effect, is controlled by 
the smaller sample dimension: the thinner the sample, the smaller the effect, while the rest of the geometry does not 
play any relevant role. The characteristic time in tr/x units is approximately equal to r/cr/i ~ 0.083 a^, when a < b, 
r/a^ ~ 0.048 b^ when b < a, and has its minimum t/<t/j, ~ 0.035 b^ for a square rod. In all cases, it stays between 3 
and 10% of the squared relevant sample size. 

Figs. I2l3i and[4]show the behavior of F, M, and r both as a function of parallel dimension b, for a = 1, and of the 
transvers dimension a, for 6=1. The straight lines are fits with the asymptotic behaviors calculated for the large a 
and b limits. 




FIG. 2: Log-log plot of the damping coefficient F in I^a units as a function of b for a — 1, and as a function of a for 6=1. 
The fits correspond to the analytic calculation. 





FIG. 4: Log-log plot of the characteristic time r in a^a units as a function of b for a = 1, and as a function of a for 6=1. 
The fits correspond to the analytic calculation. 
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VI. COMMENTS AND CONCLUSIONS 



The asymmetry in Barkhausen pulses is due to the non-instantaneous response of the field to the domain wall 
displacement. Most treatments of Barkhausen noise assume the eddy current drag on the magnetic domain wall 
to be instantaneous, which corresponds to assuming = within the magnetic domains. However, eddy currents 
take a finite time to set up after the magnetic reversal, and persist for a finite time after the corresponding wall 
displacement. This time delay causes the eddy pressure on the wall at a given time t not to be strictly proportional to 
the instantaneous velocity of the wall, but to depend on a weighted average of velocities of the wall up to time t. This 
quasi-inertial effect is responsible for the asymmetry in the pulse shape observed in Barkhausen experiments. The 
delay has a characteristic timescale t, and therefore its effect is more severe on avalanches of comparable duration. 
For very long avalanches the separation of timescales is such that the response of the field can be assumed to be 
instantaneous, therefore the asymmetry disappears, and strict universality is recovered. 

Starting from the full Maxwell equation, that includes the dynamic eddy current effects, we calculate the retarded 
pressure on the wall. The first order correction to the quasi-static (/i = 0) solution, leads to identify a damping 
coefficient, which survives to the /i — > limit and coincides with the one calculated assuming an instantaneous 
response, and an effective mass, that vanishes in the ^ —^ limit. A damping coefficient and an effective mass can 
be formally defined beyond the first order approximation, however, they both results to be frequency dependent. 
Damping coefficient, effective mass, and characteristic time depend on the sample geometry, however, it turns out 
that the only geometrical parameter that significantly affects them is the smaller sample dimension. 

The effective mass results to be negative at all frequencies. This may be understood by observing that the retarded 
pressure at time t is proportional to an average of previous velocities of the wall up to time t. When the wall is 
accelerating, the effective average velocity is smaller than the instantaneous one, and the opposite is true when the 
wall decelerates. The leftward pulse asymmetry observed in Barkhausen experiments is indeed consistent with a 
negative effective mass: the avalanche start fast, and end slowly, which is exactly the opposite of what one would 
expect from standard inertia. 

APPENDIX A: SUM OF THE SERIES Si(q) FOR a = 1 AND IN THE LIMIT a ^ oo, a ^ 

We here calculate the sum of the series 



OO 



n,fc— 1 

for some specific values of a. The value of Si for a = 1 can be obtained by writing a closed equation for Si(l) as 
follows: 



y^l 1 _y;/l l\l_y^ll y^l 1 ^^^^ 



n,A;— 1 n,fc— 1 n,/c— 1 1 

which gives 




(A3) 



so that 



Ei(l) = 7rVl28. (A4) 
To calculate Si(a) in the limit a ^ we use Eq. [iniwith z = m:/2a to write 



E°° 1 1 TT ■y^ tanh(ri7r/2ck) (A5) 



which, in the limit a ^ gives 



I]i(a^0)~^A3. (A6) 
4a 



In the opposite limit a oo, using again Eq. IA5I and the sum A„ = J2k^ ^ ~ 7r^/8 we get 

I]i(a^O)^-J^. (A7) 

APPENDIX B: SUM OF THE SERIES £2(0) IN THE LIMIT a ^ 00, a ^ 

We here calculate the sum of the series 



1 / 1 \ 



n.fc— 1 

for some specific values of a. To calculate Si (a) in the limit a ^ let us write E2 as 



11/1 1 \ 
I]2(a) = lim- y — 77^-^^ 77^^ (B2) 



,,fc=i 



and then use Eq. [Tniin both terms with z — n7r/2a and z = \J + ti: /2a respectively 



^ , , , Iv^ 1 TT / tanh(7rn/2a) tanh (tt-s/ti? + e 2a) \ ,^ , 

S2 a = hm - > [ , . (B3) 

^ ^ c^0€^n^4:ay n Vn^ + e J 



Expanding for small e and taking the limit one gets 



, , TT ■^tanh (7rn/2a) tt^ 1 — tanh^ (7rn/2a) , , 

^2(«) = ^2. -5 -4 ■ (B4) 



8a^—^ n° 16a n 

n 71 

Now we can finally take the limit a ^ to get 

n 

The limit a — > 00 can also be obtained from Eq. IB4I bv using the expansion tanh(a;) ~ x — /3 for small arguments, 
which gives 

^^'^^°°'~8a^\2an^ 3{8a)^ ny 16a2Z^U* {2a)^ J ~ 768a^ ' ^ ' 
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